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Abstract 

Without assuming that Higgs masses have the same values as other scalar masses at the 
input GUT scale, we combine constraints on the minimal supersymmetric extension of the 
Standard Model (MSSM) coming from the cold dark matter density with the limits from 
direct searches at accelerators such as LEP, indirect measurements such as 6 — » S7 decay 
and the anomalous magnetic moment of the muon. The requirement that Higgs masses- 
squared be positive at the GUT scale imposes important restrictions on the MSSM parameter 
space, as does the requirement that the LSP be neutral. We analyze the interplay of these 
constraints in the (/x, m^), (/i, mi/2), (^1/2? "^0) and (m^,tan/5) planes. These exhibit new 
features not seen in the corresponding planes in the constrained MSSM in which universality 
is extended to Higgs masses. 
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1 Introduction 



In order to avoid fine tuning to preserve the mass hierarchy mw -C mp supersymmetry 
at the TeV scale is commonly postulated. Cosmology also favours the TeV mass range, if 
the lightest supersymmetric particle (LSP) is stable, as occurs if R parity is conserved. In 
the following, we work within the minimal supersymmetric extension of the Standard Model 
(MSSM). Furthermore, in order to satisfy the strong constraints on charged or colored dark 
matter P], we require that the LSP is neutral. Over almost all of the MSSM parameter space 
the lightest neutral sparticle is a neutralinoQ i-e., a mixture of the B, W"^, Hi and H2. In 
this case, the LSP would be an excellent candidate for astrophysical dark matter 0]. 

A key uncertainty in the MSSM is the pattern of soft supersymmetry breaking, as de- 
scribed by the scalar masses mo, gaugino masses mi/2 and trilinear couplings Aq 0]. These 
presumably originate from physics at some high-energy scale, e.g., from some supergravity 
or superstring theory, and then evolve down to lower energy scale according to well-known 
renormalization-group equations. What is uncertain, however, is the extent to which uni- 
versality applies to the scalar masses ttiq for different squark, slepton and Higgs fields, the 
gaugino masses mi/2 for the SU{3), SU{2) and U{1) gauginos, and the trilinear couplings A 
corresponding to different Yukawa couplings. We do not consider here non-universal gaugino 
masses or A parameters. 

The suppression of flavour-changing neutral interactions ||^ suggests that the mo may 
be universal for different matter fields with the same quantum numbers, e.g., the different 
squark and slepton generations 0. However, there is no very good reason to postulate 
universality between, say, the spartners of left- and right-handed quarks, or between squarks 
and sleptons. In Grand Unified Theories (GUTs), there must be universality between fields 
in the same GUT multiplet, e.g., UL,dL,UR and en in a 10 of SU{5), and this would extend 
to all matter fields in a 16 of SO{10). However, there is less reason to postulate universality 
between these and the Higgs fields. Nevertheless, this extension of universality to the Higgs 
masses (UHM) is often assumed, resulting in what is commonly termed the constrained 
MSSM (CMSSM). Alternatively, there may be non-universal higgs masses (NUHM) in the 
more general MSSM. 

We and others have previously made extensive studies of the allowed parameter space 
in the CMSSM p, |10|, |ll|], incorporating experimental constraints and the requirement 
that the relic density f^x^^ ^^^^ within a range 0.1 < f^x^^ ^ ^-^ preferred by cosmology. 
There have also been many studies of the NUHM case in the MSSM ||T^, |13|, [1^, |^, in which 
the character of the LSP may change, perhaps becoming mainly a Higgsino H, rather than 
a Bino B as in the CMSSM. We think it is opportune to study the NUHM case again, 
taking into account more recent improvements in the understanding of the cosmological relic 
density Q^h!^, including x ~ ^ ^ind x ~ t [Q coannihilations in addition to x ~ x' ~ 
coannihilations [|15| and the roles of direct-channel MSSM Higgs resonances [0, as well 
as the (almost) final versions of the direct constraints imposed by LEP experiments. 



We comment below on specific cases wliere tlie sneutrino may be tlie LSP. 
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2 Experimental Constraints and the NUHM Analysis 



Important experimental constraints on the MSSM parameter space are provided by direct 
searches at LEP, such as that on the hghtest chargino x^'- '^x='= ~ 103.5 GeV [|T7|, and that 



on the selectron e: ^ 99 GeV [T^, depending shghtly on the other MSSM parameters. 
For our purposes, another important constraint is provided by the LEP lower limit on the 
Higgs mass: itlh > 114 GeV [|19[ in the Standard Model. This limit may be applicable to 



the lightest Higgs boson h in the general MSSM, although possibly in relaxed form 0. We 
recall that nih is sensitive to sparticle masses, particularly m^, via loop corrections [EO, |2T[|: 



6ml^^ ln{--^\ +... (1) 




which implies that the LEP Higgs limit constrains the MSSM parameters. 

We also impose the constraint imposed by measurements of 6 — 57 BR(i? 
X.'j) = (3.11 ± 0.42 ± 0.21) X 10"^ which agree with the Standard Model calculation 
BR{B X,7)sM = (3.29 ± 0.33) x 10"^ [g. Typically, the 6 ^ 57 constraint is more 
important for /i < 0, but it is also relevant for /i > 0, particularly when tan f3 is large. 

We also take into account the anomalous magnetic moment of the muon. The BNL 
E821 p4| experiment reported a new measurement of = |((7^ — 2) which deviates by 



1.6 standard deviations from the best Standard Model prediction (once the pseudoscalar- 
meson pole part of the light-by-light scattering contribution |^ is corrected). Currently, 
the deviation from the Standard Model value is —6 x 10^^'^ to 58 x 10~^° at the 2-cr level. 



The 2-0" limit still prefers [26] the fi > part of parameter space, but /i < is allowed so 



long as either (or both) mi/2 and tuq are large |2^. Where appropriate, the current (7^ — 2 
constraint is taken into account 0. 

In the following, we display the regions of MSSM parameter space where the supersym- 
metric relic density = flxPcriUcai falls within the following preferred range: 



0.1 < a^h' < 0.3. (2) 

The upper limit is rigorous, and assumes only that the age of the Universe exceeds 12 Gyr. 
It is also consistent with the total matter density Qm ~ 0.4, and the Hubble expansion 
rate h ~ to within about 10 % (in units of 100 km/s/Mpc). An upper limit stronger 



than (H) could be defensible |2^ , in particular because global fits to cosmological data may 
favour a lower value of flm- On the other hand, the lower limit in (^ might be weakened, in 
particular if there are other important contributions to the overall matter density. 

In the CMSSM, the values of the Higgsino mixing parameter /i and the pseudoscalar 
Higgs mass are determined by the electroweak vacuum conditions, once mi/2, mo, tan/3 
and the trilinear supersymmetry-breaking parameter Aq are fixed. This is no longer the case 

^As we discuss later, there is no such relaxation in the regions of MSSM parameter space of interest to 

us. 

■^We note that a new experimental value with significantly reduced statistical error is expected soon, but 
the theoretical interpretation will still be subject to strong-interaction uncertainties in the Standard Model 
prediction, so the impact may be muffled. 
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when the universahty assumption is relaxed for the Higgs multiplets. No longer assuming 
that the Higgs soft masses mi and m2 are set equal to rriQ at the GUT scale, we are free to 
choose /i and as surrogate parameters. Thus we have as our free parameters mo, mi/2, 
Ao, tan/?, /X and tua- The results depend somewhat on the masses of the top and bottom 
quarks that are assumed: for definiteness, we use the pole mass rrit = 175 GeV and the 
running mass mi,{mb)s^M = 4.25 GeV. 

As noted above, we include x ~ ^ = e, /i, r) and x ~ t Ell coannihilations in our 
calculation. However, since we assume for simplicity that Aq = 0, x ~ t coannihilation is 
not very important. Since fi is now a free parameter, the neutralino can be Bino-like or 
Higgsino-like depending on the ratio of /i to mi/2- When the LSP is Higgsino-like, it is often 
nearly degenerate with the second lightest neutralino x' the lightest chargino x^ [HI 



All relevant coannihilation processes for this case, have been included. In some small regions 
of the parameter space X; x'; and / can all be degenerate, and we have also taken into 
account the x' ~ ^ ^ind x^ ~ ^ coannihilations. 



3 Analysis of the (/i, m^) Plane 

We start our analysis by studying the range of possibilities in the (/x, m^) plane for various 
fixed choices of mi/2, rriQ and tan/?. Panel (a) of Fig. |l] displays the {n^rriA) plane for 
mi/2 = 300 GeV, mo = 100 GeV and tan/? = 10. The very dark (red) regions at large 
|/x| appear where the LSP is no longer the neutral x but the charged lighter stau fi, which 
is unacceptable astrophysically. There are also small 'shark's teeth' at ~ 400 GeV, 
niA ~ 300 GeV where the fi is the LSP. At large the fi is driven light primarily by the 
large mixing term in the stau mass matrix. At small particularly at small m^ when the 
mass difference m\ — m\ is small, first the fji mass is driven small, making the fi the LSP 
again. However, at even smaller |/i| the lightest neutralino gets lighter again, since m^^, ~ /z 
when /i < Ml ~ 0.4mi/2- 

The light (turquoise) shaded region in panel (a) of Fig. |l] is that for which 0.1 < Vt^h"^ < 
0.3. We see strips adjacent to the fi LSP regions, where x ~ coannihilation is 
important in suppressing the relic density to an acceptable level. The thick cosmological 
region at smaller ^ corresponds to the 'bulk' region familiar from CMSSM studies. The 
two (black) crosses indicate the position of the CMSSM points for these input parameters. 
Extending upward in uia from this region, there is another light (turquoise) shaded band 
at smaller Here, the neutralino gets more Higgsino-like and the annihilation to W^W~ 
becomes important, yielding a relic density in the allowed range 0. For smaller the relic 
density becomes too small due to the X ~ x' " coannihilations. For even smaller 
(~ 30 GeV) many channels are kinematically unavailable and we are no longer near the h 
and Z pole. As a result the relic density may again come into the cosmologically preferred 
region. However, this region is excluded by the LEP limit on the chargino mass as explained 
below. 

"'The Higgs multiplets i?i,2 couple to d,u quarks, respectively. 
'^This is similar to the focus-point region |^ in the CMSSM. 
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Figure 1: Compilations of phenomenological constraints on the MSSM with NUHM in the 
{fi,mA) plane for ta.n P = 10 and (a) nriQ = 100 GeV, mi/2 = 300 GeV, (h) rriQ = 300 GeV, 
mi/2 = 1000 GeV, assuming Aq = 0, rrit = 175 GeV and mb{mb)g^^j = 4.25 GeV. The light 
(turquoise) shading denotes the region where 0.1 < f^^/i^ < 0.3, and the (blue) solid line is 
the contour = mAf^, near which rapid direct- channel annihilation suppresses the relic 
density. The darker (green) shading shows the impact of the b ^ s'y constraint, and the 
darkest (red) shading shows where the LSP is charged. The dark (black) dashed line is the 
chargino constraint m^± > 104 GeV: lower values are not allowed. The lighter (red) 
dot-dashed line is the contour = 114 GeV calculated using FeynHiggs f^^: lower rriA 
values are not allowed. The dark (black) dot-dashed line indicates when one or another Higgs 
mass-squared becomes negative at the GUT scale: only lower and larger tua values are 
allowed. The crosses denote the values of fi and niA found in the GMSSM. 

The unshaded regions between the allowed bands have a relic density that is too high: 
O^/i^ > 0.3. However, the f coannihilation and bulk bands are connected by horizontal 
bands of acceptable relic density that are themselves separated by unshaded regions of low 
relic density, threaded by solid (blue) lines asymptoting to tua ~ 250 GeV. These lines 
correspond to cases when ~ mAfi, where direct-channel annihilation: x + x — A, is 
important, and suppresses the relic density |T6| creating 'funnel'-like regions. 

Overlaying the cosmological regions are dark (green) shaded regions excluded by 6 57. 
That for /i < is more important, as expected from previous analyses. Taking this into 
account, most of the bulk and coannihilation regions are allowed for /i > 0, but only the 
coannihilation regions and a small slice of the bulk region for n < 0. In this example, the 
'funnel' regions are largely excluded by 6 — >■ 57 for both signs of fi. 

We next explain the various contours shown in Fig. ^ The dark (black) dashed line 
in Fig. |l|(a) is the contour m^± = 104 GeV, representing the LEP kinematic limit on the 
chargino mass. The actual LEP lower limit is in fact generally somewhat smaller than this. 
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depending on the details of the MSSM parameters, but the differences would be invisible 
on this plot. We see that the chargino constraint excludes the |/i| ~ 100 GeV strip. The 
lighter (red) dot-dashed line in Fig. |l](a) is the LEP lower limit of 114 GeV on the mass of 



the lightest MSSM Higgs boson h, as calculated using the FeynHiggs programme This 



limit is relaxed in certain regions of parameter space because of a suppression of the ZZh 
coupling but, as discussed below in connection with Fig. ^, this relaxation is irrelevant for 
the conclusions presented here. For the choices of the other MSSM parameters in Fig. |l](a), 
the main effect of the Higgs constraint is essentially to exclude negative values of /x. The 
pale (pink) solid line at /i = marks the — 2 constraint. This constraint excludes the 
/i < half-plane, while allowing all of the /x > parameter space for this particular choice 
of mi/2 and uiq. 

The darker (black) dot-dashed line in Fig. |l](a) indicates where one or the other of the 
Higgs mass-squared becomes negative at the input GUT scale, specifically when either {m\ + 
/i^) < or (m2 + /i^) < 0. One could argue that larger values of and/or smaller values of 
iriA are excluded by requiring the preferred electroweak vacuum to be energetically favoured 
and not bypassed early in the evolution of the Universe. However, for a different point of 
view, see |]30[| . 

As noted above, the two (black) crosses in Fig. |l|(a) are the two pairs (/x, m^) ~ (±390, 450) 
GeV obtained in the CMSSM, assuming that the soft supersymmetry-breaking Higgs scalar 
masses are equal to the universal squark and slepton masses- squared at the GUT input scale 
(UHM). For this particular choice of the other MSSM parameters mi/2,fno, Aq and tan/3, 
these two CMSSM points both yield relic densities ^^^'^ ~ 0.21, within the range preferred 
by astrophysics and cosmology. On the other hand, the CMSSM point with /i < has 
unacceptable values for 6 — 57 decay, a^, and ruh, although the latter might conceivably 
be salvaged if the theoretical approximations in FeynHiggs happen to err in the favourable 
direction. 

The main conclusions from Fig. |l](a) are that moderate values of yU > are favoured, rriA 
cannot be small, and there is a large fraction of the remaining MSSM parameter space 
where the cosmological relic density lies within the range favoured by astrophysics and 
cosmology. This includes parts of the 'bulk' regions identified in CMSSM studies, with 
200 GeV < /i < 1000 GeV and 270 GeV < ^ 650 GeV, generalizing the CMSSM point 
for /i > 0. 

The notations used for the constraints illustrated in the other panel of Fig. |l| are the 
same, but the constraints interplay in different ways. In panel (b), we have chosen a larger 
value of mi/2- In this case, the previous region excluded by the neutral LSP constraint at 
large migrates to larger and is no longer visible in this panel. However, the 'shark's 
teeth' for moderate grow, reaching up to rriA ~ 1250 GeV. These arise when one combines 
a large value of mi/2 with a relatively small value of mo, and one may find a f or even a e 
LSP. The large value of mi/2 also keeps the rate of 6 — > 57 under control unless rriA is very 
small. The chargino constraint is similar to that in panel (a), whereas the rrih constraint 
is irrelevant due to the large value of mi/2- The — 2 constraint does not provide any 
exclusion here. Finally, we observe that the GUT mass-squared positivity constraint now 
allows larger values of ^ 1800 GeV. In this example, the two CMSSM points are at 
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(/i, myl) ~ (±1100, 1330) GeV and both have rehc densities that are too large: fi^^^ — 

The main conclusions from Fig. |I](b) are that moderate values of /i are favoured, which 
may be of either sign, but still ttla cannot be small. As in panel (a), there is a large fraction 
of the remaining MSSM parameter space where the cosmological relic density lies within the 
range favoured by astrophysics and cosmology, for both signs of fi. 

The reader may be wondering by now: how non-universal are the Higgs masses in the 
previous plots? Do they differ only slightly from universality - in which case the usual 
CMSSM results would be very unstable numerically, or do they involve violations of univer- 
sality by orders of magnitude - in which case the NUHM discussion would be unimportant 
and the usual CMSSM results would be very stable numerically? Some answers are provided 
in Fig. H, where we plot contours of rfii = sign(mf) x |mi/mo| as paler (red) curves and of 
1712 = sign(m2) X |m2/mo| as darker (black) curves. We see in panel (a), for tan/? = 10, 
mi/2 = 300, mo = 100 GeV and Aq = 0, corresponding to panel (a) of Fig. |l], that relatively 
large values are attained for both rhi and rh2, ranging up to 15 or more in modulus. How- 
ever, comparing Fig. |^(a) and Fig. ||(a), we see that such large values are not attained in 
the restricted region of the (/i, m^) plane that are allowed by the various phenomenological 
constraints discussed earlier. We have indicated by darker shading the regions of Fig. ^(a) 
which are allowed by the non-cosmological constraints, and see that |mi,2| ~ 5 in the regions 
allowed. 

Turning now to panel (b) of Fig. |, for tan/3 = 10, mi/2 = 1000, rriQ = 300 GeV and 
Aq = 0, we see values of |mi 2I ^ 5 throughout the portion of the plane displayed. Comparing 
again with the corresponding panel (b) of Fig. |I|, we see that most of the displayed range of 
11 is allowed by the experimental constraints, but only for ^ 600 GeV. In this region, we 
find |mi_2| ^ 5. 

We conclude that the answer to the questions posed earlier lie in between the extremes 
proposed. The variations in |mi,2| are by no means negligible, but neither are they excessive 
within the regions of parameter space allowed by the experimental constraints. The CMSSM 
results are not excessively unstable, but significant variations are possible for plausible ranges 
of mi,2- 

We comment finally on the magnitude of the ZZh coupling sin^(/5— a), which controls the 
observability of the lightest MSSM Higgs boson at LEP. The contours sin^(/? — a) = 0.7, 0.5 
and 0.3 are shown as closely-spaced dashed lines at the bottom of each panel in Fig. 0. In 
the bulk of the (/x, tjia) planes shown, the ZZh coupling is close to its Standard Model value, 
and the LEP lower limit on rrih is indistinguishable from the Standard Model value limit of 
114 GeV. 

4 Analysis of the (/i, M2) Plane 

We now turn to projections of the MSSM parameter space on the {fi, M2) plane, which is 
often used in the analysis of the chargino and neutralino sectors of the MSSM. In the past, 
it has been demonstrated how the cosmological constraint on the relic neutralino density 
may be satisfied in a large part of the (/i, M2) plane, for certain values of the other MSSM 
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Figure 2: Contours of the scaled Higgs masses rhi and 1712 in the {fi, m^) plane for tan /5 = 10 
and (a) rriQ = 100 GeV, mi/2 = 300 GeV, (h) niQ = 300 GeV, mi/2 = 1000 GeV, assuming 
Aq = 0, rut = 175 GeV and mb{mb)^^ = 4.25 GeV. As in Fig. the light (turquoise) 
shading denotes the region where 0.1 < ^x^"^ ^ O-^- -^^^ darker shading denotes the region 
not excluded by the other constraints. The dark (black) lines correspond to contours of 1712, 
and the lighter (red) lines to contours ofrhi. In (a) the thick solid contours correspond to 
171 = 0, the thick dashed contours to m = ±5, and the thin solid contours to m = —lb. In 
(b) the thick solid contours correspond to rh = 0, the thick dashed contours to rh = 2, the 
thin solid contours to rh = 3, and the thin dashed contours to rh = ±5. The dark (blue) 
dashed lines at very low values of indicate the contours sin^(/5 — a) = 0.7,0.5 and 0.3, 
which decrease with ttla. 
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parameters |T^, [13], |T^, ^ . Here we limit ourselves to a couple of examples that indicate how 
the cosmological region may vary. 

Panel (a) of Fig. |^ displays the (/x, M2) plane for the choices tan/3 = 10, mo = 100 GeV, 
niA = 700 GeV and Aq = 0. In this case, the region favoured by cosmology, shown by the 
light (turquoise) shading, is in the part of the plane where the LSP x is mainly a Bino, as 
preferred in the CMSSM. The LEP chargino constraint, shown as a dark (black) dashed line, 
excludes small values of fi and/or M2, where a Higgsino LSP might have constituted the dark 
matter. The rrih constraint, shown as a paler (red) dot-dashed line, excludes low values of 
M2, particularly for /i < 0, but allows substantial fractions of the cosmological regions. The 
6 ^ 57 constraint, shown in darker shading, excludes another part of the remaining allowed 
region for /i < 0, but leaves almost untouched the allowed region for /i > 0. The — 2 
constraint, shown as a pale (pink) solid line, excludes a larger region of /i < 0, but leaves 
an allowed region at higher M2. The requirement that the effective Higgs masses- squared 
be positive at the GUT scale allows a triangular region centred around /x = and extending 
up to M2 ~ 900 GeV, bounded by the dark (black) dot-dashed lines, and is compatible 
with all the other constraints in regions with both signs of fi. The cosmological region is 
bounded above by the dark (red) shaded region where the LSP is charged, close to which 
coannihilation is important in suppressing the relic density to an acceptable level. Also 
shown at large M2 is the solid (blue) line where = mA/2, but this has little effect on the 
relic density in the region allowed by the other constraints. 

Panel (b) of Fig. |^ displays the (/i, M2) plane for the choices tan P = 10, mo = 400 GeV, 
rriA = 700 GeV and Aq = 0. In this case, the cosmological region is largely complementary 
to panel (a), since both the Bino and Higgsino regions are excluded. Only narrow strips in 
the regions where the LSP is a strong mixture of gaugino and Higgsino have an acceptable 
relic density, with the exception of indentations where m^ ~ ^1^/2 and broader regions at 
large M2, where contributions from s-channel annihilation and f coannihilation are both 
important. The chargino, Higgs, h —>■ 57, Qp. — 2 and GUT positivity constraints interplay 
much as in panel (a). 

We have also examined the (/i, M2) planes for other choices of mo and m^. With a 
judicious choice of these parameters, a large fraction of the domain where /i ^ M2 and the 
LSP is gaugino-like may happen to have an acceptable relic density, though it may then be 
excluded by other constraints. However, we have not found any instance where a mainly 



Higgsino-like LSP is permitted with a sufficiently large relic density: see also [^. This is 
largely due to the LEP lower limit on m^± and the fact that large M2 is excluded at small 
by the neutral LSP requirement. 



5 Analysis of the (mi/2,mo) Plane 

This projection of the MSSM parameter space has often been used in studies of supersym- 
metric dark matter, in particular in the context of the CMSSM, where it was instrumental 
in the specification |3l| of benchmark scenarios compatible with all the experimental and 
cosmological constraints discussed earlier. Here we discuss examples in the more general 
MSSM context, which indicate some of the range of new possibilities that it opens up. 
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tan p = 10 , mo = 100 GeV, = 700 GeV tan p = 10 , nio = 400 GeV, 111^ = 700 GeV 
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Figure 3: Compilations of phenomenological constraints on the MSSM with NUHM in the 
(/i,M2) plane for tan/3 = 10 and (a) mo = 100 GeV, tua = 700 GeV^b) ttiq = 400 GeV, 
ruA = 700 GeV, again assuming Aq = 0, nit = 175 GeV and mb{mb)s^fj = 4.25 GeV. The 
pale (turquoise) shading denotes the region where 0.1 < ^^^"^ ^ ^-2' (blue) solid 

line is the contour = mA/2, near which rapid direct- channel annihilation suppresses the 
relic density. The darker (green) shading shows the impact of the b ^ s'j constraint, and the 
darkest (red) shading shows where the LSP is charged. The dark (black) dashed line is the 
chargino constraint m^± > 104 GeV: lower values of \fi\ and/or M2 are not allowed. The 
lighter (red) dot-dashed line is the contour rrih = 114 GeV calculated using FeynHiggs /|^7|/.- 
lower ruA are not allowed. The pale (pink) solid line shows the region excluded by — 2. 
The dark (black) dot-dashed triangular line indicates when one or another Higgs mass-squared 
becomes negative at the GUT scale: only lower and intermediate mi/2 are allowed. The 
(black) crosses denote the GMSSM points. 
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Figure 4: Compilations of phenomenological constraints on the MSSM with NUHM in the 
(mi/2,mo) plane for tan/5 = 10 and (a) fi = 600 GeV, ttla = 400 GeV, (h) /i = 400 GeV, 
tua = 1000 GeV, again assuming Aq = 0, nit = 175 GeV and mb{mb)s^M = 4.25 GeV. The 
pale (turquoise) shading denotes the region where 0.1 < ^x^"^ < 0-3; (^"iT-d the (blue) solid line 
is the contour = niA/'^, near which rapid direct- channel annihilation suppresses the relic 
density. The darkest (red) shading shows where the LSP is charged. The dark (black) dashed 
line is the chargino constraint m^± > 104 GeV: lower mi/2 are not allowed. The lighter (red) 
dot- dashed line is the contour = 114 GeV calculated using FeynHiggs /|17|/.- lower mi/2 
are not allowed. The light (pink) solid line is where the supersymmetric contribution to the 
muon anomalous moment is = 58 x 10^^^ : lower rriQ and mi/2 are excluded at the 2 — a 
level. The dark (black) dot-dashed lines indicates when one or another Higgs mass-squared 
becomes negative at the GUT scale: only intermediate values of mi 12 are allowed in panel 
(a), and larger values in (b). The (black) cross in panel (b) denotes the GMSSM point. 



Panel (a) of Fig. || shows the (mi/2,''^o) plane for tan/? = 10 and the particular choices 
/i = 600 GeV and tua = 400 GeV. The dark (red) shaded regions are excluded because 
the LSP is charged: the larger part resembles the similar excluded regions in the GMSSM. 
As in the GMSSM studies, there are light (turquoise) shaded strips close to these forbidden 
regions where coannihilation suppresses the relic density sufficiently to be cosmologically 
acceptable. Further away from these regions, the relic density is generally too high. The 
near-vertical dark (black) dashed and light (red) dot-dashed lines are the LEP exclusion 
contours m^± > 104 GeV and nth > 114 GeV. As in the GMSSM case, they exclude low 
values of mi/2, and hence rule out rapid relic annihilation via direct-channel h and Z° poles. 

A striking feature when mi/2 ~ 400 to 500 GeV is a 'funnel' with a double strip of 
acceptable relic density. This is due to rapid annihilation via the direct-channel A, H poles 
which occur when = niA/'i = 200 GeV, indicated by the solid (blue) line. Inside the 
shaded double strip, the funnel contains an unshaded strip where the relic density falls below 
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the range preferred by cosmology in the absence of other types of cold dark matter. The 
existence of analogous rapid-annihilation funnels has been noticed previously in the CMSSM 
context: there they were diagonal in the (mi/2,wio) plane, because the CMSSM imposes a 
link between mo and rriA that is absent in the more general MSSM discussed here. 

There is also another strip in Fig. ^(a) around mi/2 ~ 1100 GeV where the relic density 
falls again into the allowed range. In this region, the neutralino acquires enough Higgsino 
content for the relic density be in the range preferred by cosmology. For larger mi/2, the relic 
density is suppressed even more hj x ~ x' ^ coannihilation, and for mi/2 ~ 1200 GeV 
the relic density falls down below 0.1. This strip along with most of the (mi/2,'^o) plane 
shown in Fig. ^a) is actually excluded by the requirement that the Higgs scalar masses be 
positive at the input GUT scale, as indicated by the dark (black) dot-dashed lines. However, 
the rapid- annihilation funnel is still allowed by this constraint, as is a part of the x ~ ^ 
coannihilation region. The CMSSM point is located beyond the region shown here. 

Panel (b) of Fig. ^ shows the {mi/2,mQ) plane for tan/3 = 10 and the different choices 
fi = 400 GeV and rriA = 1000 GeV. In this case, the dark (red) shaded charged-LSP region 
has rather different shape, but still excludes a substantial region with large mi/2 and small 
rriQ. The striking feature of Fig. ^b) is the broad band of allowed relic density around 
"^1/2 ~ 700 GeV, where the relic density is suppressed into the preferred range by % — — 
coannihilation, the relic density falling below 0.1 for mi/2 ~ 800 GeV. In this case, the region 
allowed by the GUT stability requirement extends up to mi/2 ~ 1600 GeV, and is therefore 
satisfied throughout most of the displayed part of the {mi/2,mo) plane. 

In the small mo ~ 200 GeV and mi/2 ~ 300 GeV corner, the tau sneutrino can become 
lighter than the stau, due to the negative value of m2 — mf. Including x~ '^t coannihilation 
would shift the region preferred by cosmology for the x LSP to somewhat higher values of 
mo and mi/2- Note that, in a very narrow strip between the x- and f-LSP regions, the 
LSP is in fact u^. The direct-channel H, A rapid-annihilation region is not seen in panel 
(b), because the neutralino has already become Higgsino-like with a mass ~ 400 GeV when 
mi/2 ~ 1000 GeV, and therefore m^ is always below m^/2. For other MSSM choices, the 
rapid-annihilation region may overlap with the broad band where x~x' ^X^ coannihilation 
begins to become important. 

These two examples serve to demonstrate that the (mi/2,mo) plane may look rather 
different in the CMSSM from its appearance in the CMSSM for the same value of tan f3. In 
particular, the locations of rapid-annihilation funnels and X — x' — coannihilation regions 
are quite mo del- dependent, and the GUT stability requirement may exclude large parts of 
the (mi/2,mo) plane. 



6 Analysis of the (m^,tan/3) Plane 

This projection of the MSSM parameter space is often used in discussions of Higgs phe- 
nomenology, and a number of benchmark Higgs scenarios have been proposed As we 
now see in more detail, these do not always take fully into account other phenomenological 
constraints on the MSSM. We give examples of (m^,tan/5) planes for selected values of the 
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Figure 5: Compilations of phenomenological constraints on the MSSM with NUHM in 
the (m^, tan/3) plane for mo, ^i) = (a)(600, 800, 1000), (b)(250, 1000, 200) GeV. The 

lighter (red) dot-dashed lines are the contours rrih = 114 GeV, the solid (blue) lines show 
where ~ ^iaI'^ and the dark (black) dot-dashed line indicates when one or another Higgs 
mass-squared becomes negative at the GUT scale. The light (turquoise) shading indicates 
where 0.1 < Q^h'^ < 0.3, the darker (green) shaded regions are excluded by the b ^ s'y 
constraint. 



other MSSM parameters in Fig. ^. We first note the following general features. The LEP 
constraint on rrih excludes a region at low rriA and/or tan/5, and the b s'-f constraint 
also removes large domains of the {rriA, ta.n (3) planes, whose locations depend on the other 
MSSM parameters. The requirement that the Higgs masses-squared be positive at the GUT 
scale also excludes large domains of parameter space. However, the most striking constraint 
is that imposed by the cosmological relic density. 

We discuss first panel (a) of Fig. ^ for (mi/2,mo,/u) = (600,800,1000) GeV. This is 
similar to one of the Higgs benchmark scenarios proposed in ||3^, which has fi = 2000 GeV. 
However, we find that, for such a large value of fi, the GUT positivity constraints would be 
disobeyed all over the (myi,tan/5) plane. This reflects the fact the good renormalization- 
group running up to the GUT scale was not considered as a criterion in selecting the Higgs 
benchmark points. Reducing fi to the value 1000 GeV shown in panel (a) of Fig. |^, the GUT 
positivity constraints exclude a region at small and tan (3, which is largely excluded also 
by the b ^ s'j constraint. In this panel, there is only a very narrow range of values of 
where the relic density falls within the range favoured by cosmology. It corresponds to one of 
the 'funnels' noted previously in an analysis of the CMSSM at large tan/? 0, and is actually 
divided into two strips. The broader strip of values is just above ~ 2m^, and a 
narrower strip appears just below ttla ~ 2m^. In between, there is a narrow strip of 
values where the relic density lies below the preferred cosmological range: this strip would 
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be allowed if there is another important source of cold dark matter. Outside the double strip 
around ~ 2m^, the relic density is higher than allowed by cosmology. This exclusion 
could be evaded only by postulating that the lightest neutralino is unstable, either in an 
i?-conserving model if there is a lighter sparticle such as the gravitino, or in an i?-violating 
model. Thus most of the plane in panel (a) of Fig. |^ is excluded, even after reducing fi. 

The picture changes strikingly in panel (b) of Fig. |^, for the choices {mi/2,iTLQ, fi) = 
(250,1000,200) GeV. The contour nih = 114 GeV is shown as a lighter (red) dot-dashed 
line descending steeply when ~ 130 GeV, and then flattening out at tan/3 ~ 10. The 
vertical solid (red) line shows where rriy. ~ mA/2. The dark (green) shaded regions at 
(mA,tan/3) ^ (300 GeV, 20) and ^ (600 GeV, 45) are excluded by the 6 — S7 constraint. 
We see in this case that the relic density constraint, indicated by the light (turquoise) shading, 
is satisfied throughout a broad swathe of ^ 250 GeV for tan /? = 5 to ^ 450 GeV for 
tan/3 = 50. The requirement that the Higgs masses-squared be positive at the GUT scale 
excludes a domain of parameter space at low tan /3, as indicated by the dark dot-dashed line. 

The Higgs benchmark scenarios were not chosen with the cosmological relic density in 
mind, whereas this was taken explicitly into account in formulating the sparticle benchmark 



scenarios proposed in |^T|. It would be interesting to study NUHM Higgs benchmark sce- 
narios that respect the cosmological relic density constraint more systematically. This could 
perhaps be done by postulating a value of mi/2 that varies with tan/3. However, a detailed 
study of this point goes beyond the scope of this paper. 



7 Conclusions 

We have shown in this paper that relaxing the scalar-mass universality assumption for the 
MSSM Higgs multiplets opens up many phenomenological possibilities that were not evident 
in CMSSM studies with universal masses. We have emphasized the importance of requiring 
the LSP to be neutral and imposing the positivity of scalar masses- squared at the GUT 
scale. We find that a mainly Bino neutralino LSP is still preferred, though it may also 
be mixed with a large Higgsino component. However, we do not find MSSM parameter 
regions where the LSP is mainly a Higgsino. In addition to direct-channel A, H poles and 
X — ^ coannihilation, we have identified generic instances where X~ 'x' ~ coannihilation 
is important. We have also shown that Higgs benchmark scenarios do not respect in general 
the full range of phenomenological requirements. It is desirable to look for an agreed set of 
MSSM benchmarks that incorporate these in the studies of Higgs phenomenology, which we 
do not attempt here. 

The higher-dimensional parameter space of the MSSM with non-universal Higgs masses 
is clearly much richer than the simplified CMSSM, and we have only been able to touch on 
a few of the more striking aspects in this paper. More work is needed to digest more fully 
the impacts of the different experimental, phenomenological, theoretical and cosmological 
constraints. There are surely many more interesting features beyond those mentioned here. 
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